Astronomy & Astrophysics manuscript no. opt'thick'regions 


© ESO 2012 


February 9, 2012 





The effect of local optically thick regions in the long-wave 
emission of young circumstellar disks 

L. Ricci^'^ F. Trotta^'^^ l_ Testi^'^ A. Natta^'^, A. Isella^, and D. J. Wilner^ 
^ European Southern Observatory, Karl-Schwarzschild-Strasse 2, D-85748 Garching, Germany 

^ Division of Physics, Mathematics and Astronomy, California Institute of Technology, MC 249-17, Pasadena, CA 
91125, USA 

^ INAF - Osservatorio Astrofisico di Arcetri, Largo Fermi 5, 1-50125 Firenze, Italy 

* School of Gosmic Physics, Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin 2, Ireland 
^ Harvard-Smithsonian Genter for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 

Received XXX 2010/ Accepted YYY ZZZZ 

ABSTRACT 

Multi-wavelength observations of protoplanetary disks in the sub-millimeter continuum have measured spectral indices 
values which are significantly lower than what is found in the diffuse interstellar medium. Under the assumption that 
mm-wave emission of disks is mostly optically thin, these data have been generally interpreted as evidence for the 
presence of mm/cm-sized pebbles in the disk outer regions. 

In this work we investigate the effect of possible local optically thick regions on the mm-wave emission of protoplanetary 
disks without mm/cm-sized grains. A significant local increase of the optical depth in the disk can be caused by the 
concentration of solid particles, as predicted to result from a variety of proposed physical mechanisms. We calculate the 
filling factors and implied overdensities these optically thick regions would need to significantly affect the millimeter 
fluxes of disks, and we discuss their plausibility. 

We find that optically thick regions characterized by relatively small filling factors can reproduce the mm-data of 
young disks without requesting emission from mm/cm-sized pebbles. However, these optically thick regions require 
dust overdensities much larger than what predicted by any of the physical processes proposed in the literature to drive 
the concentration of solids. We find that only for the most massive disks it is possible and plausible to imagine that 
the presence of optically thick regions in the disk is responsible for the low measured values of the mm spectral index. 
For the majority of the disk population, optically thin emission from a population of large mm-sized grains remains 
the most plausible explanation. The results of this analysis further strengthen the scenario for which the measured low 
spectral indices of protoplanetary disks at mm wavelengths are due to the presence of large mm/cm-sized pebbles in 
the disk outer regions. 
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1. Introduction 

Planets around solar-like stars are thought to originate 
from the material contained in a circumstellar "protoplan- 
etary" disk. Observations of protoplanetary disks around 
pre-main sequence (PMS) stars at optical and infrared 
wavelengths have provided evidence for the presence of dust 
grains as large as at least a few /xm in many of these sys- 
tems. Since these grains are larger than the submicron-sized 
grains found in the interstellar medium (ISM, e.g. Mathis et 
al. I1977|) . these observational results have been interpreted 
in terms of dust grain growth from an original ISM-like dust 
population in the disk. These are the first steps of growth of 
solid particles which may potentially lead to the formation 
of planetesimals and then planetary bodies. 

In order to investigate the presence of larger grains in 
the disk, observations at longer wavelengths are needed. 
Furthermore, since the dust opacity decreases at longer 
wavelengths, whereas infrared observations are sensitive to 
emission from the disk surface layers, observations in the 
millimeter probe the denser disk midplane, where the whole 
process of planetesimal formation is supposed to occur. 



In the last two decades several authors measured rela- 
tively shallow slopes a of the Spectral Energy Distribution 
(SED; Fi, ^ ly" with a ~ 2 — 3) of protoplanetary 
disks at sub-mm and mm wavelengths (e.g. Beckwith & 
Sargent [M D Wi lner et al. [^UOTl Te sti et al. IMITl mm 
Natta et al. HQD l Wiln er et al. HOO pl And rews et al. [MI51 
Rodmann et al. 120061 Ricci et al. I2010a|) . Under the as- 
sumption of completely optically thin emission and if the 
emitting dust is warm enough to make the Rayleigh- Jeans 
approximation hold true at these wavelengths, the SED 
spectral index a is directly linked to the spectral index /3 of 
the dust opacity coefficient kJ3 through /3 = a — 2. In this 
way, the measured low values of a translate into values of 
/3 < 1 which are significantly lower than the value of 1.5 — 2 
associated to the ISM (Mathis et al. I1977p . For all the rea- 
sonable models of dust analyzed so far, the obtained values 
of f3 for young disks can be interpreted only if grains have 
grown to sizes of at least a few millimeters in the outer disk 
regions (see e.g. Draine 1 20061 Natta et al. I2007p . 



^ At these long wavelengths Hi„ is well approximated by a 
power-law. 
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The presence of mm/cm-sizcd solid particles in the outer 
disk poses serious problems to the models of dust evolu- 
tion in young circumstellar disks. Because of pressure, gas 
moves with sub-keplerian velocity, and thus it is slower 
than the solid component which is rotating at keplerian 
velocity. This generates a gas headwind felt by the solid 
particles, causing them to lose angular momentum and spi- 
ral inward in the disk. Dynamical models of solids in pro- 
toplanetary disks show that mm/cm-sized particles in the 
outer regions of disks are expected to rapidly drift toward 
the inner disk (Weidenschilling I1977p . As a consequence, 
these particles should not be present in the disk outer re- 
gions. This expectation is challenged by the observations of 
mm-slope of protoplanetary disks, if interpreted as due to 
large grains emission. Furthermore, Birnstiel et al. (I2010ap 
have recently shown that even if radial drift is completely 
halted by some physical mechanisms, the low mm-spectral 
indices 2 — 2.5) measured for the faintest disks cannot 
be reproduced by models of dust evolution in disks with 
low masses. This is because for these disks the gas densities 
in the outer regions are so low that mm/cm grains should 
be never formed there through sticking after mutual grain 
collision. 

Another a priori possible scenario for the interpretation 
of the measured low values of the mm-spectral indices is 
that a significant fraction of emission at these wavelengths 
come from optically thick regions in the disk. In this case, 
the spectral index of the SED would reflect only the spectral 
index of the Planck function, which is 2 for emission in 
the Rayleigh- Jeans regime. This value is consistent with 
that measured for a large population of young circumstellar 
disks (e.g. Rodmann et al. 120061 Ricci et al. I2010a|) . 

Under the assumption of smooth disks, the innermost 
regions of typical T Tauri disks (i? < 10 — 20 AU) can 
be dense enough to be optically thick even at long mm- 
wavelengths (e.g. Beckwith et al. [15501 Testi et al. lMJT]) . 
However, interferometric observations of disks can con- 
strain the physical structure of disks and quantify the im- 
pact of this optical-depth effect. The general conclusion of 
these studies is that disks are typically found to be rela- 
tively large, with outer radii of the order of ~ 100 AU or 
more. By modelling the mm-wave emission one finds that 
the typical contribution of the optically thick innermost 
regions to the total fluxes is nearly insignificant at wave- 
lengths > 1 mm ( see e.g. Testi et al. [20051 Wilner et al. [MI51 
Isella et ah [MM Guilloteau et al. [20TT|) . 

In the last years several different physical processes with 
the potential of concentrating particles in disks have been 
proposed as possible triggering mechanisms for the forma- 
tion of planetesimals (see Chiang Sz Youdin [2010p . These 
all lead to a local increase of the particle density in some 
regions of the disk. If the density gets high enough, these 
regions might become optically thick even at mm wave- 
lengths. Where in the disk these particle concentrations 
are supposed to occur depends on the particular physical 
mechanism driving the process, but they are not necessarily 
expected only in the innermost disk regions. 

If this scenario was viable, the spectral index of the SED 
would not carry the information on the grain sizes (through 
/3) and no constraints on that property could be derived 
from the observations. Therefore it is important to investi- 
gate the potential effect of optically thick disk regions on 
the disk total emission at millimeter wavelengths to under- 
stand whether the measured fluxes can be really used to test 



models of dust evolution in protoplanetary disks. With this 
work we want to answer the following questions: 1) could 
the observed low values of the (sub-)millimeter spectral in- 
dices of disks be explained by local concentrations of small 
particles (sizes < 1 mm) rather than by the widespread 
presence of mm/cm-sized pebbles? 2) If yes, which charac- 
teristics do they need to have? 3) Are the required local 
concentrations of small particles in disks physically plausi- 
ble? 

In Section [2] we present the sample of protoplanetary 
disks that we will use for our analysis. Section [3] outlines 
the disk and dust models considered for the analysis. The 
results are described in Section |4j whereas a discussion on 
the impact of the optically thick regions on high-angular 
resolution observations of disks in the sub-millimeter, and 
on the physical plausibility of the required structures are 
presented in Section |5j The main conclusions of this work 
are summarized in Section (5) 



2. Sample 

We describe here the sample of disks that we will con- 
sider in the following analysis. This is made of the sam- 
ples of low-mass young stellar objects (YSOs) in the 
Taurus, Ophi uchus a nd Ori on Ne bula SFRs by Ricci 
et al. (|2010a[ I2010bl I2011al I2011b[ respectively). These 
sources have been selected (1) by being low-mass Class 
II YSOs with no evidence of extended emission from a 
parental envelope, (2) by having known sub-mm/mm SED, 
through at least one detection in the 0.45 < A < 1.3 mm 
spectral range and a detection at A « 3 mm, (3) as well 
as no evidence of any stellar companion at spatial sepa- 
rations of about 10 — 400 AU that would likely affect the 
structure of the disk outer regions, to which observations 
at sub-mm/mm wavelengths are most sensitive to. 

In addition to these 47 sources, we consider here also 
the disks around BP Tau, DQ Tau, V836 Tau for which we 
obtained new CARMA observations and which satisfy the 
selection criteria described above, as detailed in Appendix 
A. 

For a detailed discussion on the completeness 
of the selected samples we refer to the Ricci et 
al. (|2010al[20T0bll2011all2011b|) papers. Here, we note that 
the second selection criterion listed above implies that the 
selected disks are relatively bright in the sub- mm and there- 
fore relatively massive. A list of the disk mass constrained 
from the sub-mm/mm SED of each disk is reported in 
Appendix B. As can been seen from Table 6, disks with 
masses as low as about 10~^ Af© in dust, or ^ 10^'^ Mq in 
gas (given a dust-to-gas mass ratio of 0.01), are included in 
our sample. 



3. Analysis 

In this Section we present the method adopted to explore 
the impact of optically thick regions in the total emission of 
a disk at sub-mm/mm wavelengths. Section [3.1 1 describes 
the disk model used in our analysis and how the optically 
thick regions are included in the disk structure. Section [3.2l 
describes the main properties of the dust models used to 
derive the dust opacities. 



Ricci et al.: Optically thick regions in protoplanetary disks 



3 



3.1. Disk structure 

We start by defining an "unperturbed" disk structure on 
the top of which we wiU add the contribution of optically 
thick disk regions. For a disk with midplanc temperature 
TmidiT) and dust mass surface density '^{r) between an in- 
ner and outer radius Rm and i?out, respectively, the sub- 
mm/mm SED can be modelled as the sum of the contri- 
bution by infinitesimal annuli with radius r (in cylindrical 
coordinate system): 



^unp _ 



d2 



B, 



.(T,nid(r))(l-e-^-W)27rrdr, (1) 



where i is the disk inclination, i.e. the angle between the 
disk axis and the perpendicular to the plane of the sk}^, d 
is the distance of the disk from the observer, and B^, is the 
Planck function. In Eq. [T]the optical depth T^{r) is defined 
as 



(2) 



where Ki, is the dust opacity coefficienlQ fSection l3.2p . and 
Sdust('') is the dust surface density. 

In this work we adopted a modified version of the two- 
layer models of passively irradiated flared disks developed 
by DuUemond et al. (j20011 which follows the schcmatization 
by Chiang & Goldreich ll997p to derive the disk structure. 
In order for these models to calculate the temperature of 
the disk midplane, a surface density profile has to be given 
as input. 

For the radial profile of the dust surface density we 
consider here the self-similar solution for a viscous disk 
(Lynden-Bell & Prinde fT974| : 
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with values of the Sq, 7, Tc parameters in the ranges 
observationally constrained by sub-arcsec angular resolu- 
tion imaging of young protoplanetary disks in the sub-mm 
(IseUa et al.HOOS Andrews et ah [^0051 I^CTTUl see Sectionil). 

After setting the dust surface density, the properties of 
the dust grains and of the central star (see below) , the tem- 
perature profile Tmidir) is derived by balancing the heating 
due to the reprocessed stellar radiation by the disk surface 
layers with cooling due to the dust thermal emission. Here 
we consider a PMS star with a mass of 0.5 Mq, bolometric 
luminosity of 0.9 Lq and effective temperature of 4000 K, 
which are typical values for the sample of low-mass PMS 
stars studied in Ricci et al. (|2010al[20T0bll2011all2011bp in 
the Taurus, Ophiuchus and Orion Nebula Cluster forming 
regions. 

On the top of the "unperturbed" disk structure just 
defined, we add the contribution of optically thick regions 



^ This parameter does not affect the sub-mm/mm significantly 
as long as the emission at these wavelengths is dominated by the 
optically thin disk regions. However this it not true anymore 
when emission from optically thick regions becomes important. 
In this paper we have adopted a value of 30°for the disk incli- 
nation of all the modelled disks. 

^ Note that here we have implicitly assumed that the dust 
opacity does not change within the disk. 



in the disk, i.e. regions where the optical depth >> 1 
even at sub-mm and mm wavelengths. By noticing that 
for a completely optically thick disk the Eq. [T] for the flux 
density reduces to 



thick 



COS I 



flout 



B^{Tn^idir))2TTrdr, 



(4) 



we modify Eq. [T] by adding at each radius r a fraction / of 
the area of the disk that is optically thick, so that the total 
flux density becomes: 



d2 



Rin 

X 2TTrdr 



B.{T^Ur)) [/ + (1-/) (l-e-^-W) 



(5) 



The "filling factor" / can thus be read as the fraction 
of area in the disk which is optically thick at all the wave- 
lengths considered in the analysis. In general / = f{r), i.e. 
the filling factor is a function of radius. If / = at all radii, 
Eq. [S] reduces to Eq. [1] i.e. no optically thick regions are 
added to the flux density of an unperturbed disk. If / = 1 
at all radii, Eq. [5] reduces to Eq. |4]for the flux density of a 
completely optically thick disk. 

In Section HTI we will treat firstly the simplest case of 
/(r) constant with radius, then we will consider for /(r) 
a family of step functions which are non-zero only within 
specified regions in the disk in Sect. 14.31 

Note that in Eq.|4]and[5]we assumed for simplicity that 
the temperature of the optically thick regions at a stellocen- 
tric radius r is equal to the one of the unperturbed disk 
calculated by the two-layer disk model at the same radius. 
Although not strictly correct, this is a very good approxi- 
mation. The temperature of the disk is roughly proportional 
to the ratio between the optical depth in the optical (where 
the disk absorbs most of the radiation energy) and the op- 
tical depth in the infrared (where the disk emits most of its 
radiation; see Dullemond et al. l200ip . Therefore, at first or- 
der, the dust temperature does not depend on the absolute 
value of the optical depth. 

Since for nearly all the disks imaged at high angular 
resolution in the sub-mm 7 is larger than (Andrews et 
al. [^niM I^HTUl Isella et al. [1003), Eq. [3] implies that the 
dust surface density typically decreases with distance from 
the central star. This means that the innermost regions 
of unperturbed disks can be dense enough to be optically 
thick at long mm wavelengths, as it was recognized since 
the very first analyses of sub-mm observations of YSOs 
(e.g. Beckwith et al. [TMI| . Testi et al. (pgUTj) quantified 
the impact of optically thick emission from inner disk re- 
gions by showing how the predicted mm-fluxes change when 
considering different parameters for the disk structure. To 
visualize their main results, they used the i^mm vs amm 
diagram, where -Fmm is the flux-density at a given mm- 
wavelength and amm is the mm-spectral index. In this pa- 
per we will extensively use this diagram. The main differ- 
ence with this and other previous works is that our atten- 
tion will mostly focus on the impact of possible optically 
thick regions throughout the disk, rather than only on the 
optically thick inner regions in the unperturbed disk struc- 
ture. A parametric study of the effect of these local optically 
thick regions on the mm-SED, and of their properties is pre- 
sented in Section 21 The physical origin of these regions is 
discussed in Section [2121 
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3.2. Dust opacity 

The only physical quantity present in Eq. [1] [5J [5] which is 
left to be described is the frequency-dependent dust opacity 
coefficient k^. This term, which at long wavelengths repre- 
sents the level of emissivity of the disk per unit dust mass, 
can be specified only when a model for the dust grain chem- 
ical composition, porosity, shape and size is considered (see 
e.g. Natta et al. lMH)) . 

In this paper we adopt the same dust grain model 
as in Birnstiel et al. (|2010ap . Ricci et al. (|2010a[ I2010ap . 
i.e. porous composite spherical grains made of astronom- 
ical silicates, carbonaceous materials and water ice (opti- 
cal constants from Weingartner & Draine 120011 Zubko et 
al. 119961 and Warren [19841 respectively). The refractive in- 
dices of the different species have been combined by using 
the Bruggeman mixing theory. We use the fractional abun- 
dances estimated by Semenov et al. (j2003p . After setting 
the chemical composition, porosity and shape of the grain, 
a grain size distribution n{a) has to be specified in order 
to determine the dust opacity coefficient of the dust 
model. We adopt a power-law grain size distribution with 
power-law exponent q 

n(a) oc a~'^ (6) 

truncated between the minimum and maximum grain sizes 
Omin and Omax, respectively. Since in this work we want to 
investigate the impact of local optically thick regions in the 
outer disk on the slope of the mm-SED without invoking the 
presence of large mm-sized grains, we consider here a size 
distribution which does not contain mm-grains, i.e. with 
Omax ~ 0.1 mm. Note that, for any given grain chemical 
composition, as long as Cmax ^0.1 mm the dust opacity 
coefficient Ki, at mm wavelengths stays unchanged and does 
not even depend on q (see Fig. 3 in Ricci et al. I2010ap . 
Therefore, the dust model adopted here refers to a more 
general dust population in which only grains smaller than 
^0.1 mm are present. 

4. Results 

In this section we present the results of our analysis. In 
particular, in Sect. 14.11 and 14.31 we compare the mm-flux 
densities as derived by our disk models for different val- 
ues of parameters related to the disk str ucture, with th ose 
measured by Ricci et al. (I2010al I2010bl I2011al I2011bp for 
49 disks in the Taurus, Ophiuchus and Orion star forming 
regions. 

4.1. Case of a constant f : the -Fimm-ai-3mm diagram 

Wc start by considering a /(r)-function which is constant 
for stellocentric radii lower than 300 At0 and equal to 
outside 300 AU. Each panel in Figure[l]shows the prediction 
of our disk models for different values of the filling factor 
/ and of the dust mass for a given couple of the (7,rc) pa- 
rameters for the disk unperturbed structure^. The values of 
the parameters considered in this work for the unperturbed 
disk, i.e. 7 = 0, 1 and Tc = 20, 200 AU (see equation[31), lay 

* For all the disks considered in this work, more than 80% of 
the unperturbed disk mass is contained within 300 AU. 

^ The value of the normalization factor in equation [31 i.e. Eo, 
is set after fixing the (unperturbed) dust mass in the disk. 



at the limits of the ranges for these parameters as recently 
constrained by a high-angular resolution survey of about 
30 protoplanetary disks in the sub-mm (Isella et al. 120091 
Andrews et al. I^OIM I^IJTU)) . 

In each panel the points with / = represent the emis- 
sion of the unperturbed disk, i.e. without any addition of 
optically thick regions. The dependence of the flux density 
at 1 mm on the dust mass is due to the fact that the bulk 
of the material reside in the outer disk regions where the 
surface density is relatively low and the emission optically 
thin. In terms of the spectral index ai^smm the points with 
/ = lay at relatively large values of about 3 — 3.5. This 
is due to the adopted dust model made of small grains 
fSection l3.2p . for which the value of the spectral slope (3 of 
the dust opacity coefficient is about 1.5. The slight decrease 
of Q!i-3mm with increasing dust mass is given by the fact 
that for disks that are massive enough the innermost disk 
regions (r < 10 — 20 AU) are so dense that the emission 
from these regions becomes optically thick in the sub-mm, 
and this has the effect of making the SED shallower. 

Note however that the points with / = do not get 
into the area of the i^imm-ai-3mm diagram which contains 
the bulk of the observational data. In particular the models 
presented here with / = overpredict the observed spectral 
index ai^snim- 

The effect of increasing the filling factor / of optically 
thick regions has always the same kind of signature in the 
(^iinm,ai-3mm)-plane: the absolute flux in the (sub-)mm 
increases because more and more optical depth is added 
into the system, whereas the spectral index tends to de- 
crease and approach the value of about 2, as expected for 
completely optically thick emission in the Rayleigh- Jeans 
regim^EI. 

The main result shown in Figure [T] is that, for most 
of the unperturbed disk structures, relatively low values 
of the filling factor (/ < 0.05) are required to explain the 
bulk of the data. The most significant variation of the mm- 
fluxes with the properties of the unperturbed disk comes 
when the characteristic radius changes from 20 to 200 AU 
(see left and right panels, respectively). With Tc = 20 AU 
the surface density is more concentrated in the inner disk 
than in the Tc = 200 AU case. As a consequence of this 
the impact of the optically thick inner regions to the total 
emission is more evident (see e.g. the low ai-smm-value 
of the model with / = and a dust mass of 10 ~^ Mq 
for the unperturbed disk structure). However, the fluxes of 
the disks with / > do not dramatically depend on the 
structure of the unperturbed disk. This is due to the fact 
that in most of the cases already for / w 0.05 the total disk 
emission becomes dominated by the added optically thick 
regions. 

What varies strongly with the unperturbed disk struc- 
ture at a given / is the mass which has to be included in 
these added regions in order to make them optically thick. 
This is discussed in the next section. 



Note that in some cases the disk models predict values for the 
spectral index ai_3mm which are lower than 2. This is because 
the Rayleigh- Jeans regime is not completely reached at these 
wavelengths in the outer and cooler regions of the disks. 
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Fig. 1. Flux density at 1mm vs spectral index between 1 and 3mm. In each panel red, blue and green points represent 
observational data of disks in Class II YSOs in Taurus (Ricci et al. l2010a1 and this work), Ophiuchus (Ricci et al. l2010b|) . 
and Orion Nebula Cluster (Ricci et al. I2011al I2011bp . respectively. The 1 mm- fluxes in the Orion Nebula Cluster have 
been multiplied by a factor of (420 pc/140 pc)^ to compensate for the different distances estimated for the Orion (~ 
420 pc, Menten et al. [MJ7| and Taurus/ Ophiuchus regions (- 140 pc, Bertout et al. HM^ Wilking et aL [^nU5| . Red 
empty circles are the three Class II YSOs observed with CARMA and discussed in the Appendix A. Black symbols show 
the model predictions for different unperturbed dust mass: triangles for 10~^ M©, stars for 10~^ Mq, and squares for 
10~^ Mq. Solid lines connect disk models with the same unperturbed dust mass, whereas dashed lines connect models 
with the same constant filling factor / for the optically thick regions, with values indicated in the bottom right side of 
each panel. In each panel a given couple of the (7, rc)-parameters has been assumed for the unperturbed disk structure, 
with values indicated in the top left corner. 



4.2. Case of a constant f: required masses in optically thick 
regions 

So far we have considered the emission of our modelled 
disks without considering how much mass has to be added 
to the disk in order to make a fraction of the disk surface 
optically thick even at millimeter wavelengths. A simple 
way to do this is to calculate the surface density which has 
to be present in the added regions to give them an optical 
depth of about 1 at the longest wavelength considered in 
this work, i.e. 3 mm, and then integrate over the disk area 
covered by the filling factor /. This represents actually a 
lower limit for the mass which has to be present in the 
added regions to make them optically thick: if more mass 
is put onto those regions this does not have a significant 
effect onto the SED. 



Therefore, the surface density of dust in the optically 
thick regions S^yg^. is given by the condition 

T3mm ~ ^iust'^Smm ~ 1 Sjugt ~ 1/ KSmm, (7) 

which implies that Sj^^gj does not depend on the stellocen- 
tric radius. As a consequence of this the total amount of 
dust in these regions 

<.t - / ^Lt^^rdr « f^Lt^Rlu (8) 

i.e. the total dust mass in the optically regions depends 
quadratically on the largest radius in which these regions 
are present in the disk, namely 300 AU in this simulation. 
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For the dust considered in Section [321 '<3mm ~ 0.45 cm^/g 
and M^^^^ « / x 0.07 A/©. This means that even in the 
case of the model with the largest unperturbed mass in 
dust (10^^ -^^o) and with the lowest value of the filling 
factor / considered here (0.01), the optically thick regions 
contain at least 70% of the mass in the unperturbed disk. 
This ratio then increases linearly with / and decreases with 
increasing the dust mass of the unperturbed disk. This ar- 
gument clearly shows that in order to have optically thick 
regions of the kind discussed so far in the outer disk a very 
strong redistribution of dust particles has to occur in the 
disk (see Section [5]). 

As shown in Eq. |51 the dust mass in the optically thick 
regions depends quadratically on the disk outer radius. An 
idea to reduce the amount of dust required in these optically 
thick regions is therefore to reduce the area in the disk 
where they are located. This is the topic of the next section. 

4.3. Case of f{r) as a step function: the Fi,nm-tti-3mm 
diagram 

We consider here the result of optically thick regions lo- 
calized in smaller areas of the disk. In particular, wc dis- 
cuss four cases in which /(r) is a step function with val- 
ues greater than between 10 and 20 AU, 30 and 40 AU, 
50 and 60 AU, 80 and 90 AU, respectively. We consider 
annuli with a width of 10 AU because the physical mech- 
anisms proposed so far to concentrate particles locally in 
the disk typically act on these length scales or smaller (see 
Section [5^ . The different central radii chosen for the an- 
nuli allow us to investigate how the location of the optically 
thick regions in the disk can affect its total sub-mm emis- 
sion. 

The four panels in Fig. [5] show the model predictions on 
the i^imm — ai-3mm diagram for disks with such localized 
optically thick regions and an unperturbed disk structure 
with 7 = 0.5 and Tc = 60 AU, which are approximately the 
median values for disks imaged at high-angular resolution. 

The optically thick regions can have a signficant impact 
onto the global mm-SED even if they are concentrated in 
rings with the relatively small width of 10 AU. The effect 
is stronger for the disk with lower masses because of the 
higher constrast in optical depth between the unperturbed 
disk and the optically thick regions. For the same reason, 
i.e. the contrast with the unperturbed disk structure, and 
for the fact that the area of an annulus of a given width 
scales linearly with the central radius of the annulus, at a 
given filling factor / inside the ring, the effect of the op- 
tically thick regions is the largest in the ring which is the 
furthest from the star, i.e. between 80 and 90 AU. Note 
that to reproduce the bulk of the data, larger fractions / 
are needed than in the /(r) = const case. This is because 
in the case of optically thick emission, apart for the tem- 
perature, it is the surface area of the emitting material that 
determines the amount of its emission. Therefore, if these 
regions are distributed over a smaller area of the disk, they 
need to occupy a larger fraction of that area, which is what 
we find with our analysis. This argument justifies our choice 
of dealing with very simple structures for the optically thick 
regions: even if optically thick regions in real disks would 
likely have more complex structures than modelled here, 
our analysis is meaningful in terms of the fractional area 
covered by such regions throughout the disk. 



4.4. Case of fir) as a step function: required overdensities 

As done in Section 14.21 for the case of a constant / within 
300 AU, we analizc here the dust mass which has to be 
present in the optically thick regions. This can be calculated 
by using the first equality in Eq.[5]with the inner and outer 
radii of the ring for and i?out, respectively. 

For the case of the ring between 10 and 20 AU the 
total dust mass in the optically thick regions is M^^^^ k, 
/ X 2.4 • 10-4 Mq, and rises to w / x 5.5 • lO"'' Mq, 
« / X 8.7 • 10-"* Mq and « / X 1.3 • lO^^ M© when the 
ring is moved outward to 50-60 AU and 80-90 AU, respec- 
tively, because of the increased area of the ring itself. To 
explain the left end of the Fimm-ai-3mm diagram (Fig. [2]) 
filling factors / > 0.1 on the top of the lowest mass disk 
(Mdust.unpert = 10"'^ Mq) are needed. However, the mass in 
the optically thick regions is larger than the one in the un- 
perturbed structure for all the rings considered here. This 
would require an extremely strong concentration of par- 
ticles in those regions which are not seen in the results 
of the numerical simulations run so far (see discussion in 
Section 15. 2|) . For the disk with A^dust,unpcrt = 10"^ Mq 
a significant decrease of the spectral index ai-smm, i-C- 
down to about 2.5 and lower is obtained only for / > 0.3 
(Fig. (2). These filhng factors require dust masses in the 
optically thick regions as large as at least 70% of the un- 
perturbed disk mass. Finally, the added optically thick 
regions in the most massive disk considered here, with 
-A/dust, unpcrt = 10~^ Mq, coutain relatively low dust mass 
as compared with the unperturbed disk mass. For example, 
in the case of / = 0.5 the ratio among the former and lat- 
ter masses is about 12%, 28%, 44%, 65% for the optically 
thick regions inside rings with radii of 10-20 AU, 30-40 AU, 
50-60 AU and 80-90 AU, respectively. 

Tables[T]|4]report, for models with different unperturbed 
disk structures in terms of total dust mass, 7 and rc, the 
required overdensities in the optically thick regions. These 
are defined as the ratio between the dust surface density 
^dust optically thick regions and the surface den- 

sity in the unperturbed disk structure evaluated at the 
center of the ring. Since S^^^^ is always the same (i.e. 
l/^Smm ~ 2.2 cm^/g), the required overdensities depend 
only on the value that the dust surface density of the un- 
perturbed disk assumes in the ring, and therefore on the 
7, T-c) parameters which define such surface 
density. Note that the values of these overdensities do not 
depend on the adopted value for the dust opacity Ksmm (or 
equivalently on the adopted dust model). This is because 
both the surface density in the optically thick regions, and 
the surface density in the optically thin unperturbed disk 
structure at a given flux depend on the dust opacity as 
'^Smm- The lowest values for the overdensities are found 
for the most massive disk when the overdensity region is 
located closest to the star. The low overdensity values of 
order of unity indicate that this disk is massive enough 
to have the unperturbed structure at radii of ^10-20 AU 
marginally optically thick by itself, i.e. without the addition 
of any artificial optically thick region. 
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Fig. 2. As in Fig. [U but with optically thick regions localized in annuli of 10 AU-width, with inner and outer radii 
indicated in the top left corner of each panel. The unperturbed structure of the disk is characterized by 7 = 0.5 and 
rr = 60 AU. 
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Table 1. required dust overdensities in the added optieally thick regions inside a ring between 10 and 20 AU from the 
central star. 





7 = 1 


7 = 0.5 


7 = 


rc = 20 AU 


r-c = 200 AU 


rc = 20 AU 


rc = 200 AU 


rc = 20 AU 


Tc = 200 AU 


-A'-fdust,unpcrturb = 10 Mq 


100 


520 


70 


1200 


55 


3200 


^'^dust.unpcrturb ~ 10 Mq 


10 


52 


7 


120 


5.5 


320 


A/dust, unporturb ~ 10 Mq 


1 


5.2 


0.7 


12 


0.55 


32 



Table 2. required dust overdensities in the added optically thick regions inside a ring between 30 and 40 AU from the 
central star. 





7 = 1 


7 = 0.5 


7 = 


rc = 20 AU 


rc = 200 AU 


rc = 20 AU 


rc = 200 AU 


rc = 20 AU 


rc = 200 AU 


Afdust.unpcrturb = 10 " Mq 


630 


1300 


530 


2000 


680 


3300 


A/dust,unpcrturb = 10 Mq 


63 


130 


53 


200 


68 


330 


A^dust.unpcrturb = 10 Mq 


6.3 


13 


5.3 


20 


6.8 


33 



Table 3. required dust overdensities in the added optically thick regions inside a ring between 50 and 60 AU from the 
central star. 





7 = 1 


7 = 0.5 


7 = 


rc = 20 AU 


rc = 200 AU 


rc = 20 AU 


rc = 200 AU 


rc = 20 AU 


rc = 200 AU 


-A-^dust.unpcrturb = 10 Mq 


2700 


2300 


5700 


2700 


62000 


3400 


A/dust, unpcrturb = 10 Mq 


270 


230 


570 


270 


6200 


340 


A/dust, unpcrturb = 10 Mq 


27 


23 


57 


27 


620 


34 



5. Discussion 

5.1. Impact of the optically thick regions on high-angular 
resolution observations of disks in the sub-mm 

So far we have discussed the possible effect of local opti- 
cally thick regions in the integrated fluxes of young disks 
in the millimeter. Since the required overdensities needed 
to make these regions optically thick are typically very 
large (especially for the low-mass disks, see Tables [T]|3]) 
one could expect these regions to be easily detectable with 
high-angular resolution imaging through (sub-)millimeter 
interferometry. However, the contrast between the surface 
brightnesses inside and outside the optically thick regions 
is not the only factor determining the possible observabil- 
ity of these structures. In fact, if the bright optically thick 
regions were uniformly distributed throughout the disk and 
with characteristic length scales much smaller than the an- 
gular resolution of the observations, they would be smeared 
out by the convolution with the resolution element. 

Since all the high-angular resolution observations of 
disks conducted so far in the (sub-)mm have revealed a 
disk structure which is essentially homogeneous, the re- 
sults of existent observations provide an upper limit to the 
characteristic length scales of the invoked optically thick re- 
gions. The highest angular resolutions achieved so far in the 
sub-mm are about 0.2-0.4 arcsec, corresponding to physi- 
cal scales of ^ 30 — 60 AU at the distances of nearby star 
forming regions (Isella et al. 120101 Andrews et al. 120111 
Guilloteau et al. I201ip . Note that this is larger only by 
a factor of a few than the width of 20 AU considered in 
Section 14.31 This means that if optically thick regions were 
concentrated in those 20 AU-wide annuli, a radial bump 
in the surface brightness map could have been marginally 
detected, although not spatially resolved, by these observa- 
tions. 

In order to probe non-homogeneous structures at 
smaller scales higher angular resolution is needed. The 



Atacama Large Millimeter/Submillimeter Array (ALMA) 
will allow to do that down to scales of a few AU only. For 
example, ALMA can detect spiral density waves in nearby 
massive disks (Cossins et al. I2010p . As detailed in the next 
subsection in these regions the concentration of particles 
might be efficient enough to make the dust emission from 
spiral arms optically thick. 

If non-homogeneous regions will be detected by future 
observations with very high angular resolution, a prediction 
of the models presented in this work is that the millimeter 
spectral index measured in the bright regions should be 
equal to the spectral index of the Plank function in that 
region, which would be 2 in the case of Raylcigh- Jeans 
emission. To determine more precise predictions from these 
models would require further modelling of the disk struc- 
ture with a proper treatment of the radiative transport in 
the optically thick regions. This goes well beyond the scope 
of this work. The future results coming from observations 
with ALMA and the EVLA arrays have the potential to 
guide future developments of disk models along these lines. 

Recently, Guilloteau et al. (|201ip presented an in- 
vestigation of the radial variation of the dust properties 
through a dual- frequency (1.3 and 3 mm) survey of disks 
in the Taurus- Auriga region with the Plateau de Bure 
Interferometer. In most disks that are resolved at the two 
wavelengths, they find a radial dependence of the millimeter 
spectral index which they interpret as due to a variation of 
the spectral index of the dust opacity coefficient (3 through- 
out the disk. In particular the /3-index appears to typically 
increase from low values (^ 0-0.5) within ^ 50 — 100 AU 
from the central star to about 1.5-2 at the disk outermost 
regions. In the scenario presented in this work, the same 
data could be potentially explained in terms of a different 
efficiency of the physical mechanism(s) concentrating parti- 
cles in different regions of the disk: in the inner disk, filling 
factors of optically thick regions / ^ a few x 0.1 could 
explain the lower spectral indices (see Fig. [5]), thus mim- 
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Table 4. required dust overdensities in the added optieally thick regions inside a ring between 80 and 90 AU from the 
central star. 





7 = 1 


7 = 0.5 


7 = 


rc = 20 AU 


r-c = 200 AU 


rc = 20 AU 


rc = 200 AU 


rc = 20 AU 


Tc = 200 AU 


-A'-fdust,unpcrturb = 10 Mq 


19000 


4100 


370000 


3800 


2.4-10'' 


3700 


^'^dust.unpcrturb ~ 10 Mq 


1900 


410 


37000 


380 


2.4-10'' 


370 


A/dust, unporturb ~ 10 Mq 


190 


41 


3700 


38 


2.4-10' 


37 



icking the effect of a lower /3; in the less dense outermost 
regions the same mechanisms might be not efficient enough 
to concentrate particles and that could explain the large 
spectral indices which are typical of optically thin emission 
from small dust grains (sizes < 0.1 mm). 

From the existent data only it is not possible to rule out 
the scenario described in this paper. We also note here that, 
if present, high levels of overdensities like those invoked by 
this analysis can probably occur only in the disk midplane. 
Because of the lower optical depths at longer wavelengths, 
only observations in the millimeter (or longer wavelengths) 
can probe directly the existence of these overdense regions. 
Observations in the near and mid-infrared probe only the 
surface layers of the disk, and they would be therefore insen- 
sitive to concentrations of particles in the midplane. Also 
observations in the far-IR (^ 100 ^m) are still significantly 
affected by optically thick emission from the unperturbed 
disk structure, at least for typical disks around T Tauri 
stars with disk masses larger than a few Jupiter masses (see 
e.g. the analysis by Andrews & Williams I2005p . This would 
make the analysis presented in this paper much more model 
dependent at IR wavelengths than at longer mm-waves be- 
cause of the impact of disk structure (e.g. flaring) and the 
particular choice for the treatment of the radiative transfer 
in the disk. Note that this might be not the case, especially 
in the far-IR, for disks with lower masses e.g. around young 
brown dwarfs. At the same time, depending on the physi- 
cal mechanism driving this local concentration, the vertical 
distribution of particles in the disk, and therefore also in 
the surface layers, might be somewhat affected. In order 
to accurately determine this effect one would need global 
dynamical models of dust particles in a gaseous disk cou- 
pled with radiative transfer calculations to predict the disk 
emission in the IR. 

In the next subsection we discuss the possible mecha- 
nisms proposed in the literature to concentrate particles in 
the disk. In particular we aim to understand how plausi- 
ble are the overdensities required to explain optically thick 
regions considered in the last Sections. To do this we base 
the discussion on our current theoretical knowledge of the 
potential processes driving the concentration of solids in 
young disks. 



5.2. Physical plausibility of the required overdensities 

To understand whether the overdense optically thick re- 
gions discussed in this paper are physically plausible in 
real protoplanetary disks, one has to compare the re- 
quired filling factors and overdensities of regions with small 
< 0.1 mm-sized grains discussed in the last section with 
the outcome of the numerical simulations which investigate 
different mechanisms leading to the concentration of solid 
particles. These mechanisms are often invoked in the lit- 
erature for their potential of halting locally the otherwise 



fast radial drift of solids (Brauer et al. I2007P and of forming 
planctesimals in the disk (see a recent review by Chiang & 
Youdin l^rmUl . 

A promising mechanism to concentrate solid parti- 
cles is through the development of streaming instabili- 
ties, in which the gas and solid components are mutually 
coupled by drag forces in a turbulent disk. Johansen & 
Youdin (|2007|) showed that overdensities of particles even 
larger than 1000 can be formed. Looking at Tables [TJU] this 
indicates that these overdense regions would probably be 
optically thick also in the millimeter for most of the disk 
models considered here. However, in terms of the impact 
to the global disk SED, the filling factor / of these re- 
gions is probably not large enough to be really significant: 
the fraction of particles in overdense regions at level of 
the order of 100 or above is only about 1% or less (see 
Fig. 11 in Johansen & Youdin [2007|) . This means that the 
filling factor / of disk area occupied by such regions is 
much lower than that, as the mean distance between par- 
ticles decreases with increasing density. Furthermore, all 
these highly-overdense regions are obtained for solid par- 
ticles which are only marginally coupled to the gas, with 
values for the Stokes parameteiQ, St^l. Grains with sizes 
< 0.1 mm, as those considered in this analysis, are instead 
very well coupled to the gas (St<<l). In the case of stream- 
ing instabilities Johansen & Youdin (I2007P showed that al- 
ready for particles with St~0.1, the level of overdensity is 
decreased by a factor of 10. For these reasons, we argue 
that, although streaming instabilities may potentially form 
overdense regions which are optically thick at millimeter 
wavelengths, they would not reach the required filling fac- 
tor of the disk surface to explain the observations. 

Another possible mechanism which has been proposed 
to trap solids involves the presence of large anticyclonic vor- 
tices in the disk. These structures can be the result of baro- 
clinic instability (Klahr & Bodenheimer 12005)) . the Rossby 
wave instability (Lovelace et al. 119991 Regaly et al. I201ip . 
or magneto-rotational instability (Fromang et al. I2005p . 
Recent 2D numerical simulations of circumstellar disks have 
shown that, although dust overdensities larger than 2-3 or- 
der of magnitudes can be obtained in vortices covering a 
significant fraction of the disk surface area, the solids that 
are trapped are larger than about 10 cm, whereas much 
smaller sub-mm dust grains are not significantly affected 
by such structures (Lyra et al. I2009p . 

Long-lived axisymmetric pressure bumps have been ob- 
tained in simulations of magneto-rotational turbulent disks 
(Johansen et al. I2009p . These pressure bumps have the po- 
tential of trapping solids which are marginally coupled to 



^ The Stokes parameter is an adimensional parametrization 
for the particle size in a gaseous disk. In the Epstein regime, 
which is relevant for the conditions in the disks treated in this 
work, St= 7raps/(2Sgas), where a and ps are the particle size 
and density, respectively. 
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the gas, like nim-cm sized pebbles, but the obtained over- 
densities in gas are not large enough (<10 — 20%) to redis- 
tribute more gas-couplcd small grains at the levels required 
for the optically thick regime (Tables [T] . 

Finally, some concentration of dust particles can oc- 
cur in disks undergoing gravitational instabilites. In these 
disks, the non-linear evolution of the instabilities lead to 
the formation of spiral waves with local overdensities in 
the gas component as high as about 100 (Rice et al. 120041 
Boss l20l0|) . Since sub- mm sized particles are well coupled to 
the gas, the same level of overdensity is expected for small 
grains as well. Considering that spiral waves are character- 
ized by very extended structures, the results presented in 
this work show that the overdense regions in spiral waves 
can be optically thick even at mm- wavelengths and can even 
dominate the emission of a young disk at these wavelengths. 
This means that for gravitational unstable disks the mea- 
sured low values of the mm-spectral index (ai_3mm ^ 3.0) 
can be potentially explained by the optically thick emis- 
sion of small grains from overdense regions in spiral waves. 
However, in order for disks to devclope these instabilities, 
they need to be rather massive. Recent numerical simula- 
tions of Boss (|2010p have shown that more than about 0.04 
Mq inside 20 AU are needed around a 1 A/0 young star, 
or about 0.02 Mq around a 0.5 Mq young star as consid- 
ered in this work. By assuming a standard ISM-likc value 
of 100 for the gas-to-dust ratio, this lower-limit corresponds 
to 2 ■ 1O~^M0 of dust mass inside 20 AU from the central 
protostar. Therefore, only the mm-fluxes of disks in the 
bright tail of the i^imm — c«i-3mm diagram can be explained 
by accumulation of unprocessed small grains via gravita- 
tional instabilities. High sensitivity and angular resolution 
imaging with ALMA in the sub-mm will soon constrain 
the occurrence of these instabilities in real disks (Cossins 
et aL l^OTU)) . 

6. Summary 

In this work wc investigated the effect of possible local op- 
tically thick regions on the mm-wave emission of proto- 
planetary disks. The main goal of this exercise was to see 
what kind of structures are needed to explain the low val- 
ues of spectral indices measured for protoplanetary disks in 
nearby star forming regions, without invoking the presence 
of mm-sized particles in the disk outer regions. 

To quantify this potential effect we considered a simple 
disk model in which the disk is filled with regions with op- 
tically thick emission even at long millimeter wavelengths. 
In our calculations these regions occupy a fractional area 
of the disk which is left as a free parameter of the models, 
and wc derive the predicted mm-fluxes for a range of differ- 
ent values for the parameters defining the structure of the 
unperturbed disk. 

Interestingly, our analysis shows that relatively small 
filling factors for these optically thick regions can reproduce 
the measured fiuxes of young disks in the (sub-)millimeter. 
In most cases, if the filling factor is kept constant through- 
out the disk, values lower than a few percent are sufficient 
to match the observed data. If the regions are instead local- 
ized in annuli with a width of 10 AU, the required filling 
factors increase to a few x 10%. 

A significant local increase of the optical depth in the 
disk can be caused by the concentration of solid particles, 
as predicted by different physical mechanisms potentially 



acting in the disk. We discuss the physical plausibility of 
the required optically thick regions. In particular we inves- 
tigate whether the required overdensities arc reproduced in 
the simulations of the physical processes proposed in the 
literature to drive the concentration of particles. 

The main conclusion of this work is that for the vast 
majority of disks no physical processes proposed so far are 
capable to reproduce the measured low mm-spectral indices 
via a concentration of small (<< 0.1 mm) particles in opti- 
cally thick regions. According to our current knowledge of 
these processes, only for the brightest disks, likely the most 
massive ones, the total mm-wave emission can be strongly 
affected by optically thick spiral arms and innermost dense 
regions. The results of this analysis further strengthen the 
scenario for which the measured low spectral indices of pro- 
toplanetary disks at long wavelengths are due to the pres- 
ence of large mm/cm-sized pebbles in the disk outer regions. 
Observations of disks at (sub-) millimeter and centimeter 
wavelengths can therefore be used to constrain the models 
of the early phases of planetesimal formation. 

Appendix A: New CARMA observations at 3mm 

We describe here new mm-wave observations of nine 
YSOs in the Taurus- Auriga star forming region (SFR). 
The sources were selected for being relatively bright at 
sub-mm/mm wavelengths, namely with a flux density at 
0.85 mm greater than 100 mJy and/or a 1.3 mm-flux den- 
sity greater than 30 mJy. This criterion was chosen to have 
high chances to detect our sources at about 3 mm with 
the Combined Array for Research in Millimeter Astronomy 
(CARMA). The observed YSOs are listed in Table 1. 

The dust thermal emission toward our sample of nine 
young disk systems in Taurus- Auriga was observed with 
CARMA between 2010 March 2 and April 6. The array 
was in C configuration which provides baselines between 
30 and 350 m. Observations were carried out at a central 
frequency of 102.5 GHz (2.92 mm), with a total bandwidth 
of 4 GHz. The raw visibilities for each night were calibrated 
and edited using the MIRIAD software package. Amplitude 
and phase calibration were performed through observations 
of the QSOs 0336-f323 (for 04113+2758, V892 Tau, FN 
Tau), 3C 111 (for IC 2087, AB Aur, BP Tau, V836 Tau), 
0449-1-113 (for DQ Tau). Passband calibration was obtained 
by observing the QSO 3C 84. Mars and Uranus were used 
to set the absolute flux scale. The uncertainty on CARMA 
flux calibration is typically estimated to be 15% and is due 
to uncertainties in the planetary models and in the cor- 
rection for atmospheric effects and hardware instabilities 
(see Ricci et al. 2011 for an analysis of the repeatability 
of 3mm-flux measurements with CARMA on observations 
performed in the same time period of the ones presented 
here). Maps of the dust continuum emission were obtained 
by standard Fourier inversion adopting natural weighting, 
and photometry was obtained in the image plane. The re- 
sulting FWHM of the synthesized beam is about 2". 

Among the nine observed YSOs, wc considered the three 
sources (BP Tau, DQ Tau, V836 Tau) which satisfy the 
selection criteria adopted in this work and listed in Sect.jH 
We instead did not include the other observed sources for 
the following reasons: 

— AB Aur and V892 Tau, since they are Herbig Ae/Be 
stars, and therefore with a central PMS star more mas- 
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Object 


RA 


DEC 




rms 


Comments 


(1) 


(J2000) 
(2) 


(J2000) 
(3) 


(mJy) 
(4) 


(mJy) 
(5) 


(6) 


AB Aur 


04:14:47.8 


26:48:11.1 


10.0 


0.9 


Her big Ae/Be 


VSQ9 Tan 




97-55-fl5 5 


54.9 


1.0 


rTprni cf A p /Rp 

J-J-Cl Ulti JiC / JJC 


BP Tau 


04:30:44.3 


26:01:23.9 


10.1 


0.6 




DQ Tau 


04:38:28.6 


26:10:49.7 


14.9 


0.6 




FN Tau 


04:47:06.2 


16:58:43.0 


4.1 


0.8 


extended envelope emission 


IC 2087" 


04:39:55.8 


25:45:02.0 


< 2.7 


0.9 


Class I 


V836 Tau 


04:14:13.5 


28:12:48.8 


6.6 


0.5 




MHO 1 


04:14:17.0 


28:10:56.5 


50.5 


1.1 


no sub-mm info 


MHO 2 


04:14:47.8 


26:48:11.1 


30.2 


1.1 


no sub-mm info 



Notes. For this undetected source the (RA.DEC) coordinates are from 2MASS (Cutri et al. l2003p . 
Column (6) reports the reason why the source has not been considered in the analysis (see Appendix A) . 



sive than the sample of T Tauri stars considered in this 
paper; 

— FN Tau shows evidence of spatially resolved neb- 
ulosity in near-infrared scattering light (Kudo et 
al. I2008P , likely due to a leftover protostellar envelope. 
Furthermore, Momose et al. (j2010p showed how disk 
models fail to reproduce at the same time SMA 0.88 mm 
observations at high-angular resolution and single dish 
data at 1.3 mm, confirming that an extended "halo" 
component is present in the FN Tau system. 

— IC 2087, which is an embedded Class I YSO (Luhman et 
al. I2010p . and therefore, a significant contribution from 
an extended envelope to the sub- mm/mm emission is 
expected to occur; 

— MHO 1 and MHO 2, which are part of a binary system 
with angular separation of 4.0" (or about 560 AU at the 
Taurus distance of 140 pc, Duchene et al. I2004p and no 
intcrfcromctric observations which can separate the two 
components have been carried out so far in the sub-mm. 

As for the three selected sources, according to Luhman 
et al. (|2010p . and references therein, BP Tau is a K7- 
spectral type PMS star, DQ Tau is a double lined spectro- 
scopic binary which consists of two PMS stars with similar 
specctral type in the range of K7 to Ml, and V836 Tau is a 
K7. Andrews & Williams ([2005)) reported fluxes for BP Tau 
of 130±7 mJy and 47±0.7 mJy at about 0.85 and 1.3 mm, 
respectively, whereas for V836 Tau Fo.ssmm = 74 ± 3 mJy, 
^i.3mm = 37 ± 6 mJy. The other selected source, DQ Tau, 
is a circumbinary disk which is known to undergo recur- 
ring millimeter flares due to star-star magnetic reconnec- 
tion events (see Salter et al. 120101 and references therein) . 
By measuring the light curve in the millimeter Salter et 
al. (|2010p estimated a quiescent level of the emission at 
about 1.3 mm of 97 mJy, and of 17 mJy at 2.7 mm, which is 
roughly consistent with the measured flux of 14. 9 ±0.6 mJy 
at 3.2 mm presented in this work. By combining these 
data with the new ones obtained at 3.2 mm with CARMA 
(Table for the three selected sources, we obtained the 
mm-spcctral index shown as red empty circles in Fig. I1I2I 



Appendix B: disk properties of the selected sample 

In Table 6 wc report some of the properties of the 50 disks 
considered in this analysis. These properties were derived 
by fitting sub-mm data using two-layer disk models (for 
more details see the Ricci et al. I2010al [20T0al I2010al I2010al 
papers). 
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Table 6. Columns (2) and (3) report the source flux density at 1.0 mm and spectral index between 1 and 3 mm, 
respectively, from the best fit two-layer disk models adopted in Ricci et al. (|2010all2010b[ 1201 la[ 12011b)) . The uncertainties 
of fluxes in the millimeter range are typically dominated by the uncertainty on the absolute flux scale of the observations. 
Values of these uncertainties depend on the particular instrument and conditions during the observations, but have typical 
values around 10 — 20%. Contrary to Figures 1 and 2, the 1 mm-fluxcs reported here do not consider any correction for 
the different distances of the three regions. Column (4) lists the product between dust mass and dust opacity at 1 mm 
from the same models. Columns (5) and (6) report the inferred dust mass for values of 2.5 and 3, respectively, for the 
power-law slope q of the grain size numberdensity (see e.g. Ricci et al. I2010al for more details). Blank fields in Columns 
(4), (5), (6) are for cases in which no good fit could be obtained. 
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